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Abstract 

Let represent a two— dimensional isotropic elastic body. We consider the prob- 
lem of determining the body force F whose form ip{t){fi{x), f2{x)) with ip be given 
inexactly. The problem is nonlinear and ill-posed. Using the Fourier transform, the 
methods of Tikhonov's regularization and truncated integration, we construct a regu- 
larized solution from the data given inexactly and derive the explicitly error estimate. 
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1. Introduction 

Let n = (0, 1) X (0, 1) represent a two— dimensional isotropic elastic body. For eacli 
X := (xi,X2) G ^i, we denote hy u = {ui{x,t), U2{x,t)) the displacement, where Uj is the 
displacement in the Xj— direction, for all j G {1, 2}. As known, u satisfies the Lame system 
(see, e.g., [HIS]) 

= fiAu + {X + n)V {div{u)) + F 

where F := (Fi, F2) is the body force, div{u) = V • u = dui/dxi + du2/dx2^ and A, are 
Lame constants. We shall assume that the boundary of the elastic body is clamped and the 
initial conditions are given. 

In this paper, we shall consider the problem of determining the body force F. The 
problem is a kind of inverse source problems. The inverse source problems are investigated 
in many aspects such as the uniqueness, the stability and the regularization. There are 
many papers devoted to the uniqueness and the stability problem. In [7J, Isakov disscused 
the problem of finding a pair of functions (ti, /) satisfying 

cutt - Au = f 
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where / is independent of t. He proved that using some preassumptions on /, from the 
final overdetermination 

u{x, T) = h{x) 

, we get the uniqueness of {u, /). 

As shown in [9J, the body force (in the form (l){t)f{x)) wih be defined uniquely from an 
observation of surface stress (the lateral overdetermination) given on a suitable boundary 
of 17 X (0,T). In the paper, the authors also gave an abstract formula of reconstruction. 

Another inverse source problem is one of finding the heat source F[x, t, u) satisfying 

ut- Au = F. 

The problem was considered intensively in the last century. The problem with the final 
overdetermination was studied by Tikhonov in 1935 (see [8|). He proved the uniqueness of 
problem with prescribed lateral and final data. In the last three decades, the problem is 
considered by many authors (see O [H \TT\ \T2\ [T3| [T^ ) . Although we have many works on the 
uniqueness and the stability of inverse source problems, the literature on the regularization 
problem is quite scarce. Very recently, in [3l H] , the authors considered the regularization 
problem under both the lateral and the final overdetermination. The ideas of using the 
Fourier transform and truncated integration in the two papers are used in the present 
paper. We also consider the regularization problem under the final data and prescribed 
surface stress. 

To get the lateral overdetermination, some mechanical arguments are in order. Let 
o"i,(j2,r be the stresses (see [HE]) defined by 

T - ( ^ + dU2 \ 

^\dx2 dxij 
du ' 

(jj = Xdiv(u) + 2fi—^, j £ {1, 2} 
dxj 

We shall assume that the surface stress is given on the boundary of the body, i.e., 

/ (71 r\ / nA ^ / XA 
\r (J2) \n2) \X2) 

where X = {Xi,X2) is given on dVt, and n = (ni,n2) is the outward unit normal vector of 
dVL. 

As discussed, our problem is severely ill-posed. Hence, to simplify the problem, a 
preassumption on the form of the body force is needed. We shall use the separable form 
force as in [9] 

(Fi(x, t), F2(x, t)) = ^mfiix), f2{x)) 

where 99 is given inexactly. The form is issued from an approximated model for elastic wave 
generated from a point dislocation source (see, e.g., |9l[r0]). But, since 93 is inexact, our 
problem is nonlinear. Morever, the problem is still ill-posed because the measured data is 
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not only inexact but also non-smooth. 

Precisely, we consider the problem of identifying a pair of functions {u, /) satisfying the 
system: 

= fiAuj + (A + f^)g^,div{u) + Vj G {1, 2} (1) 

for (x, t) X (0, T), where n, A are real constants satisfying /x > and A + 2/x > 0. 
Since the boundary of the elastic body is clamped, the displacement u = {ui,U2) satisfies 
the boundary condition 

{ui{x,t),U2{x,t)) = (0,0), xedn (2) 
In addition, the initial and final displacement are given in 

' (ui(x,0),U2(x,0)) = {uoi{x),uo2{x)) 
(^(x,0),^(a;,0)) = {ul,{x),u*02{x)) (3) 

_ {ui{x,T),U2{x,T)) = {uTlix),UT2{x)) 

Finally, the surface stress is given on dO, 



(4) 



{nicTi + n2T = Xi 
7X2(72 + niT = X2 

We shall assume that the data of the system (1) — (4) 

I=i^,X,uo,u*o,UT) G {L\0,T),iL\0,T,L\dmML'mMLHn)f,iL\n)f) 

are given inexactly since they are results of experimental measurements. The system (1) — (4) 
usually has no solution; moreover, even if the solution exists, it does not depend continously 
on the given data. Hence, a regularization is in order. Denoting by lex the exact data, which 
are probably unknown, corresponding to an exact solution {uex, fex) of the system (1) — (4) 
, from the inexact data approximating I^x, we shall construct a regularized solution 
approximating fex ■ 

In fact, using the Fourier transform, we shall reduce our problem to finding the solu- 
tions of the binomial equations whose binomial term is an entire function (see Lemma 1). 
In this case, the problem is unstable in the neighborhood of zeros of the entire function. 
The zeroes can be seen as singular values. Using the method of Tikhonov's regularization 
and truncated integration, we shall eliminate the singular values to regularize our problem. 
Error estimates are given. 

The remainder of the paper is divided into two sections. In Section 2, we shall set some 
notations and state our main results. In Section 3, we give the proofs of the results. 

2. Notations and main results 
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We recall that = (0, 1) x (0, 1). We always assume that the data / = ((^, X, uo,ut, u^) 
belong to 

{L\o,T),{L\o,T,L\dn))f,{LHn)f,{LHn)f,{LHn)f) 

For all C = (6, 6), C = (Ci, C2) G we set ^ C = 6C1 + 6C2 and |C| = VFI- 
We first have the following lemma. 

Lemma 1. If u e {C^{[0,T]; L'^{n)) n L'^{0,T; H^{n)))^ , f G {L'^{n))^ satisfy (1) - (4) 
corresponding the data I, then for all a = {ai,a2) G R^\{0}, we have 

2D{I). J fj{x). cos(a • x)dx = gj{I), Vj G {1, 2} 
ci 

where 

D{I) = Di{I).D2{I),gj{I) = -^iajD2{I)ho + Diil)hj) 

with 

T T 

Di{I) = J (p{T - t)sm{^/xT2jl\a\t)dt,D2{I) = J v{T-t) sui{^\a\t)dt 


{a ■ Uq). cos(q; ■ x)dx 

+-\/A + 2fi.\a\. J {a ■ ut)- cos(q; • x)dx 
n 

- VA + 2^.|a|. cos( VA + 2n\a\T). J (a ■ uq). cos (a • x)dx 

n 

T 

- J J sin(A/A + 2/i|a|(r - t)){a ■ X). cos(a • x)dLodt 
an 

h,(i) = - si„(ViJ|«|r). / - . „;)).cos(„ . x)<fc 

+V^.H./(|a|V.-a.(«-UT)).cos(a..)dx 

— -y/^-lo:!. cos(-y//I|Q!|T). J (\afuoj — aj{a ■ uo)).cos{a ■ x)dx 

n 

T 

-J J sm{^\a\{T - t))i\a\'^Xj - aj{a ■ X)).cos{a ■ x)dudt,'ij £ {1,2}. 



an 
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From Lemma 1, we consider the function 

T T 

D{I) = J ip{T - t) sin(v^A + 2/x \a\ t)dt. J ip{T - t) sm{^\a\t)dt 



The problem is unstable in the neighborhood of zeros of this function. However, from the 
properties of analytic function, we can show that if ^ then this function differ from 
for almost every where in R"^. Furthermore, using the idea of Theorem 4 in [5j, we get the 
following lemma. 

Lemma 2. Let T,q be positive constants, ipQ £ L^(0,T)\{0} and D{ipq,t) : ^ R 

T 

D{ipo,T){a) = / (fo{t) sin{^/T\a\t)dt 



Then D{ipQ, r) 7^ for a.e a £ R? . Moreover, if we put 

q ln(£-i) 
^^-9ir-ln(ln(e-i))' > ° 

then the Lehesgue measure of the set 

Bi^o,T,e) = {a £ B{0,Re),\D{ipo,T){a)\ < e'^} 

is less than R~^ for e > small enough, where B{d,Rs) is the open ball in R? . 

Lemma 1 and Lemma 2 imply immediately the uniqueness result. 

Theorem 1. Letu.u* e {C^{[Q,T]-L'^{Vt))r\L'^{Q,T-H'^{Vt))f, fj* e {L'^{n)f. If{u,f), 
{u*,f*) satisfy (1) — (4) corresponding the same data I, and 92 ^ 0, then 

(n,/) = K,r) 

Let ( 

Uex7 fex) be the exact solution of the system (1) — (4) corresponding the exact data 
lex = iVex, Xex,UQ^ ,Uq^^ ,u^). Notice that, if we assume 

Uex e {C^{[0,T];L\n))nL\[0,T];H\n)))\fex G {L\n)f,^ex G L^(o,r)\{o} (5) 
then for all j G {1,2}, 

9jihx) 



F{fjex){ci) = 2 j fjex{x)cos{a ■ x)dx 



D{hx) 
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for a.e a G i?^, where Qj, D are defined by Lemma 1, fjex : ^ R is defined by fjex{x) 
x{^)fjex{x) + x{~^)fjex{—x), and F is the Fourier transform in R^. 
From approximate data = (if, X,uq,Uq,ut) satisfying 

II II ^ 1 1 V \r€.X 1 1 1 1 1 1 J, 



< e, \\uTj - WTillii(f^) < Vi e {1,2} 



(6) 



we construct a regularized solution /g = (/ie,/2e) whose Fourier transform is 
^(/je)(a) = x{B{Q^Re)). ^/^'}'^}!'X ya € i?2\{0} 



where 



1 J 75 



9er'Zn(/n(e-i)) 



(7) 



We have two regularization results. 



Theorem 2. Let {uex,fex) be the exact solution of the system (1) — (4) corresponding the 
exact data I^x, o-nd hold. Then from the given data satisfying (0), we can construct 
a regularized solution fi, G (C(i7))^ such that 

lim Wfj, - fjex\\L2(n) = 0, Vj G {1, 2} 



//we assume, in addition, that f^x G {H^ {^))'^ ,then 

II/,. - /,e.|li.(^) < 63er (66 ||/,ex.||^.(^) + (27r)-2^ ln(ln(6-i)) 



■ ln(e 



-1^ 



VjG{l,2} 



for e > small enough. 



Theorem 3. Let {uex,fex) be the exact solution of the system (1) — (4) corresponding the 
exact data lex, o,nd ^ hold. We assume, in addition, that 



fjex{x). cos(a • x)dx 



da < oo, Vj G {1,2} 



Then from the given data Is satisfying (0), we can construct a regularized solution f^ G 
(C(J7))^, which coincides the one in Theorem 2, such that 



fjex\ 



0, VjG{l,2} 



3. Proofs of the results 



Proof of Lemma 1 
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Proof. Let a = (01,02) G and G = cos(a • x). Notice that the j— th equation of the 
system ([1]) can rewrite 

Getting the inner product (in L'^{Q)) of the equation and G and using the condition ([2]), 
for {j, k} = {1,2}, we get 

A J UjG = j {njaj + nkT)Gduj - J (^j^dx - J -^^dx + ip{t) J fjGdx 



dn n n 



(8) 



= J XjGduj — fj^lal"^ J UjGdx — {X + fM)aj J {a ■ u)Gdx + (p{t) J fjGdx 
an n an 

Multiplying ([8]) by Uj, then getting the sum for j = 1, 2, we obtain 

J {a-u)Gdx = j {a ■ X)GdLO - {\ + 2n)\a\'^ J {a ■ u)Gdx + (p{t) J {a ■ f)Gdx (9) 
n an n n 

Multiplying ([8]) by |ap and multiplying ([9]) by —aj, then getting the sum of them, we have 

J ^|a|^ Uj — Oj. (a • u)^ Gdx = J (^\a\^ Xj — aj. {a • X)^ Gdx 

^ (10) 
—fj, |a|^ / (^\af Uj — aj. (a • n)^ Gdx + ip{t) / ^|a|^ fj — aj. (a ■ /)^ Gdx 

n n 

We consider ([9]) and (jlOp as the differential equations whose form 

y" + v^y = h{t) (11) 

where r/ is a real constant and y(0), y'{0), y{T) are given. Getting the inner product (in 
L2(0,r)) of ([II]) and sin(r/(T - t)), we have 

T 

- y'{0)sin{7]T) + rjy{T) - rjy{0)cos{rjT) = J - t) sm{r]t)dt (12) 



Applying (fT2]) to ([9]) with ry = y^(A + 2//)|q| and y = f {a ■ u).Gdx, we get 



Diil). ia-f).Gdx = ho{I) (13) 



f2 
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where Di(I),/io(/) are defined by Lemma 1. 

Similarly, applying ([12]) to ([10]) with 7] = ^/JI\a\ and y = f {\a\'^Uj — aj.{a ■ u)).Gdx, we get 

n 

D2{I). j (|q|2/,- - aj{a ■ f)).Gdx = hj{I), Vj G {1, 2} (14) 
a 

where D2{I),hj{I) are defined by Lemma 1. 

Multiplying (fT3]) by ajD2{I) and multiplying (fTi|) by Di{I), then getting the sum of 
them, we obtain the result of Lemma 1. □ 



and (/. : C ^ C 



Proof of Lemma 2 

Proo/. Put ^0 : ^ ^ -R 

( Mt) te{o,T) 
ifoi-t) te{-T,o) 

,0 t^{-T,T) 

T 

e-'^'^o{t)dt= j e-'^'^oit)dt 

-oo -T 

Then (/> is an entire function and D{ipo,T)(a) = i(j){^/T\a\) . Because (fo ^ 0, its Fourier 
transform (in R) does not coincide 0. Therefore, there exists zq £ R such that |i;A(-Zo)| = 
Ci > 0. Thus 0^0. Since is an entire function, its zeros set is either finite or countable. 
Consequently, D{ipo,T){a) ^ for a.e a G E?. 

To estimate the measure of B{ipQ, r, e), we shall use the following result (see Theorem 4 
of $11.3 in [6J). 

Lemma 3. Let f{z) be a function analytic in the disk {z : \z\ < 2eR}, |/(0)| = 1, and let 
rj he an arbitrary small positive number. Then the estimate 

ln|/(z)| > -ln(^).ln(Mj(2e/?)) 



is valid everywhere in the disk {z : \z\ < R} except a set of disks {Cj) with sum of radii 
'^rj < r]R. Where Mf{r) = max|/(2;)|. 

\z\=r 

Returning Lemma 2, we put <j)i : C ^ C 



(p{z + zo) 
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Then is an entire function, (t)i{0) = 1, and for all z E C, \z\ < 2eR, 



Ci\Mz)\ = 



1 

I 

f-T 



1 

Ke'^^.J \lpo{t)\dt = e'^^yo\\L^io,T) 

-T 



For £ > small enough, applying Lemma 3 with R= ^R^ and ij = g^^, we get 



3lnRe + ln(-^) + ln(15e^) 



.eTRs + In 



/ llyo|lLi(0,T) \ 
I ^1 ) 



> -^T.RelnR, > -gln(e-i) -ln(Ci) =ln(|^) 

for all l^l < except a set of disks {B{zj,rj)}j^j with sum of radii Y^ri < ijR = 

Consequently, for £ > small enough, we have \zq\ < ^R^ and |i?!>(2;)| = Ci. \(f)i{z — zo)\ > 
£5 for all \z\ < Re except the set U B{zj + zo,rj) . Hence, B{ipQ,T,e) is contained in the 

set U Bj, where 

Bj = {ae B{0, R,), I - % | < r,} 

with yj = Re{zj + zq). 

If Vj > ^prRe + Tj then Bj = 0. If yj < rj then Bj C 5(0, ^), so m{Bj) < If 
fj < Vj < y/rRe + Tj then 



B, C 5(0, ^4^)\i?(0, 



hence 



m{Bj) < ^(^J + ^-^^^ _ ^(^J ~ ^j)^ = ^^^J^i < 47r(V?ii:e + rj)rj 



Thus we get 



(u( MVrRe + rj)rj ^ 47rr| 
m(5(v9, r, £)) < 2^ ^ + 2^ ^ 



r — ' r 



1 



+ _.( v )^ < — 



for £ > small enough. The proof of Lemma 2 is completed. 



□ 



Proof of theorem 1 
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Proof. Put w = u — u* and v = f — f* then {w, v) satisfies (1) — (4) corresponding the data 

/= (99, (0,0), (0,0), (0,0), (0,0)) 

Let Vj : ^ R he defined by Vj{x) = x(0)fj(x) + x{—^)i^j{—x). Lemma 1 implies that, 
for all j G {1, 2}, for all a G R'^\{0}, we get 

D{I).F{vj){a) = 2D{I). j Vj{x) cos(a • x)dx = gj{I) = 

n 

Applying Lemma 2 with ipo{t) = ip{T — t), we get D{I) ^ for a.e a ^ R^. Therefore, 
Fivj) = 0, and it implies that Vj = 0. Thus v = (0,0). Hence, w satisfies that 



— - = n/\w + (A + //) V {div{w)) 
Getting the inner product (in (L^(r2))^) of p3]) and dw/dt, we have 



(15) 



1 A V 



dt 



2 

L2(n) 



i=i 



A + /X d 



2 



— \\diviw)\\L2^^^ 



Integrating this equality in (0, t), we get 

2 



E 



9t 



2 2 

+ ^ ^ II V?i;j ||^2m) + (A + //) ||dit;(u;) ||^2(f^) = 



(16) 



for all t G (0, T). Using the condition ([2]), we have 



\\div{w) Wli^^-^ = 

2 







dwj 
dxj 


L2(n) V 



dwi dw2 
dxi dx2 



= E 



+ 2 



L2(n) 



dwi dw2 
dx2 dxi 



dwi 



dx-2 



+ 



dxi 



'illL2(C) 



Since /i > and A + 2/i > 0, the above inequality implies that 

2 

f^Yl ll^'^illL2(n) + (^ + /^) N^^^("')llL2(n) ^ 
i=i 

From (fT6]) . we obtain dw/dt = (0,0). Since ?x)(x,0) = (0,0), the proof is completed. □ 
To prove two main regularization results, we state and prove some preliminary lemmas. 
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Lemma 4. Let (uex, lex) be the exact solution of (1) — (4) corresponding the exact data lex 
satisfying @, and the given data satisfying (0). Using notations of Q, we put 

G,{Ie)=x{B{0,Re)). 9j{Ie)D{Ie) 
' 5e + {D{h)f 

Then for all j G {1,2}, we have Gj{Is) G L^(i?^)nL^(i?^); moreover, there exists a constant 
Co depend only on lex such that for all e G {0,e~'^), 



Gjils) - F{fjex) < xiBiO, Re))CoRee 



l-6q 
2 



+2x{Be)\\f,ex\\L2^n)+XiR\B{0,R,)) F{f,ex, 

where B, = {a G B{0,R,), \D{Iex){a)\ < e^i} . 

Proof. First, we show that there exists a constant C2 > depend only on Iqx such, tlicit for 
aU e G (0,e-^), r > ro = q/{9T), j G {1,2}, 

\\D{Iex)ho^(^R2) < C2, \\D{Q - D{Iex)\\L^^R2) < ^2^ 
\\9jiIex)\\L'^(^B{0,r)) ^ ^'s?', \\gjile) " 9j ilex)\\ lo. ^B(0,r)) ^ ^'a^e 

Recall that Di(I), D2{I), hQ{I), hj{I) are defined by Lemma 1. For all a G i?^ we have 

\DkiIex)\ < \\(Pex\\Li{0,T) AB>k{Is) - Dkilex)\ < W^s - ^ex\\L^O,T) 

for all k G {1,2}. Hence, \D{Iex)\ < ||V5ez|lii(o,T) ^nd 

\DiIe) - D{hx)\ = \Di{h) - Di{Iex)\ ■ \D2{Ie)\ + \Dl{Iex)\ • |i?2(4) " D2{Iex)\ 
^ S-dlVexll Li{0,T) + ^) + 1 1 Vex II Li(0,T) — II V^ex || Li(0,T) 

+ e-').e 

A straightforward calculation show that, for all a G i?(0, r)\{0}, we have 

\oijhQ{Iex)\ <C:ir\af' ,\aj{hQ{Ie) - hQ{Iex))\ <Czr\a\^e, 
\h,{Iex)\ <C:ir\a\\\hj{h)-hj{Iex)\<C^r\a\'e 
for all j G {1,2}, where C3 is a positive constant depending only on lex' Therefore, 

\9Ahx)\ < l^2(/e.)| + \Dr{Iex)\ < 2C, ||9.e.|Li(0T) ^ 

a a 



and 



\a\ \a\ 

|q!| |a| 

< Gsre. (||v'ez|lii(o,T) + ^) + Csr.e + Csre. (||9'ez|lii(o,r) + ^) + Csr.e 

< 2C73 (^||(^ex||ii(o,T) + e"' + l) re 
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Returning Lemma 4, for all j e {1,2}, we get Gj{Ie) G L^{R'^) fl Lp^^R^) because the 
support of Gj{Ie) is contained in B{0,R^) and Gj{I^) G U^{R^). Moreover, 

gj{I,)D{I,) gj{h^)D{h, 



+xiB{0,R,)) 



Qj (lex) D{Iex) 9j [h. 



+ X{R\B{0,R,)). Fifjex) 



5e + {D{Ieo:)f D{h 
We shall estimate each of the terms of the right-hand side. Wc have 



9j{Ie)D{Ie) 93 [lex) Djle 



< 



6e\93 {Is) D{Ie) - 9j {lex) D{Iex)\ 

(4 + {D{h)f) (5, + {D{h.)f) 



, \D{Ie)\ . |g(Jex)| ■ \9j jle) Djle.) ' 9j {h.) D{Ie)\ 

(5, + {D{I,)f) (4 + {D{h,)f) 

^ \9j{Ie)D{Ie)-9j{Ie.)D{hx)\ \gj jh) Djlex) ' 9j jh.) DjQl 
- 6, ^ 6, 

If £ G (0, e"^) then R^ > tq, so for all a G B{0, R^) we get 

\9j{I,)D{I,) - gj{Iea,)D{Iea,)\ 

< \9j{h) - gj{Ie,)\ . \D{h)\ + \gj{Ie,)\ . \D{h) - D{Ie.)\ 

< C2Ree.{C2 +e) + C2R,e < {C2 + l)^R,e 

and similarly, 

\gj{I,)D{Iea,) - gj{h,)D{h)\ < (C2 + l)'i?.£ 
Consequently, for all e G (0, e~^), we can estimate the first term 



xiB{0,R,)) 



gj{I,)D{I,) gj{Ie^)D{Ie: 



5, + {D{h)f 5, + {D{h,)f 
Considering the second term, we have 

9j{hx)D{Iex) 9j{Iex) 



6, + {D{h,)f D{Iea^) 



<X{BiO,Re)). 
Se\9j{Iex)\ 



2{C2 + VfReE 



(5, + {D{h,)f).\D{h,)\ 



We always have 

4 \9j {Iex)\ 



< 



9j {h 



D{Ie 



(5, + {D{h:c)f)-\D{Ie.)\ 

Furthermore, if a G 5(0, Re)\Bs then 

\9j {hx)\ 



fjex{x)cos{a ■x)dx 



<2||/, 



< 



\9j {I 



(Se + {Dihx))^).\D{Iex)\ \D{Iex)\ 

12 



ex; I ^ 5£C2Re 



3 — 



-6g 



Therefore, for all £ G (0, e we can estimate the second term 

9j ihx) D{Iex) gjih 



X{BiO,Re)) 



5, + {D{I,,)f D{Ie,) 
Thus, for all £ G (0, e~^), we have 



< 2x{B,) ||/,e.|L.(n) +X(S(0,i2e))^^ 



Gj{Ie) - Fifjex) <x{B{0,Re)) 



2{C2 + lfRee , 5eC2R, 



+ 



£6g 



+2x{Be) \\fje.\\L^^^)+X{R\B{Q,R,)) F{fje 



(iq+l 



Choosing Se = £ 2 and Co = 2{C2 + 1) + C2, we complete the proof. 



□ 



It is obvious that, for all j G {1,2}, by Lebesgue's dominated convergence theorem, 



xiR'\BiO,R,)) 



jex) 



converges to in L [R ) when e ^ 0. However, to get an explicitly 



estimate for it, some a-priori information about f^x must be assume. 

Lemma 5. Let a e R, Q be an measurable subset of R^ (n> 1), and w G L^{Q) n L?'{Q). 
Then 



RP 



w{x) sin(a + akX}.)dx 



Q 



k=l 



da = 2"-i7r" \\wf 



Proof. We first prove in the case a = 0. Put w : R^ ^ R 

wix) = x{Q)w{x) - x{-Q)w{-x) 
Then w G ^^(i?") n ^^(i?") and 



/n 
w{x) sin(^ a^Xkjdx 



where F„ is the Fourier transform in i?". Using Paserval equality, we get 

2 



/n 
w{x) sin(^^ akXk)d: 

r, k=l 



L2(i?") 



(2vrr ,,^,,2 



cyn-l n ||„,,||2 

^ ^ II^IIl2(Q) 



Similarly, we also have 



Q 



w{x) cos(^^ akXk)doi 
k=l 



da = 2"-i7r" \\w"^ 
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Now, we notice that 

/n 
w{x) sm(a + akXkjda. 



(cos(a))^ 



w{x) sin(^^ akXk)d3. 



+(sin(a))^ 



fe=i 



/n 
w{x) cos(^ akXk)d. 

Q 



where 



/n „ n 

w{x)sm.(^^akXk)dx. / w{x) cos('y^ akXk)d. 
Q *^=^ n '==^ 

Since v{—a) = —v{a) for all a G i?", we get J v{a)da = 0. Thus 



w{x) sin(a + akXk)da 



k=l 



da 



= (cos(a))^2"-V" \\w\\l,^Q^ + (sin(a))^2"-V" \\w\\1,^q^ = 2^-\- \\w\\1,^q) 
The proof is completed. 

Using Lemma 5, we have the following result. 
Lemma 6. Let w e H^{n) and r > 7r/(2V2). Then 

2 

w{x) cos (a • x)dx 



Proof. Since 



R2\B(0,r) 



B?\B{0,r) 



w(a;) cos(a • x)dx 



, ^ 72V27r ,. ,,2 
da < \\w\\ fj-i 



da < J J w{x) cos(a • x)da 

^^^\aj\>r/V2 ^ 



da 



, the proof will be completed if we show that, for all j G {1,2}, 



\aj\>r/V2 



J J w{x) cos(q; • x)da 



24V27r /„ „2 
< — - — lFllL2(n) + 2 



dw 



dxj 



□ 
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We will prove for the case j = 1, and the other cases are similar. We have 

1 



J w{x) cos(q: • x)dx = J 



w{x) 



sm{a ■ x) 



xi=l 



dxo 



J xi=0 



sm[a ■ x , 



dw 



dxi ai 



dx 



so 



cos (a • x)dx 



3 



a 



w{l, X2) sin(Q;i + a2X2)dx2 



+ - 



■w{0, X2) sui{a2X2)dx2 



+ 



at 



dw 
dx\ 



. sin (a • x)dx 



Therefore, 



\ai\>r/^/2 



J w{x) cos{a ■ x)dx da < J J ^^{x).sm{a ■ x)dx 



da 



00 1 



+ J -^dai. J J w{l,X2)sm{ai + a2X2)dx2 



\ai\>r/V2 -°° 

00 1 



dan 



+ 



J -^dai. J J ■w{0,X2)sm{a2X2)dx2 



\ai\>r/V2 



-00 



da2 



127r^ 



dw 



dxi 



^ 6\/27r 
+ 

L2(n) 



1^(1' ■) IIl2(o,i) + lk(0, .) |li2(o,i) 



Noting that 



1 1 
w{l,X2) = J -^{xiw{x))dxi = J ^w{x) + xi^^{x)^ dxi 



, we get 



\w{l,X2)\'^ < J ( 2|w;(a;)p + 2 



dw 
dxi 



(x) 



dx\ 



Hence, 



1 

J \w{l,X2)f dx2 <2\\w\\l2^ii)+2 



dw 



dxi 



L2(n) 
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Similarly, 



1 i i 

J \wiO,X2)\^dX2 = j j 











d_ 

dxi 



((1 — xi)w{x)) dxi 



dxi 



1 1 / 
< j j {2\w{x)\^ + 2 

0^ 

Thus, we have 



— (x) 
dxi ^"^^ 



dxidx2 = 2 ||'w||2^2(f^) + 2 



dw 



dxi 



L2(n) 



|ai|>r/V2 
+ 



cos(a • x)dx 



da < 



dw 



dxi 



+ 



n 

24\/27r /„ ,„c^ ^ 



L2(n) 



9xi 



< 



24^/27r 



k(i,-)lli2(f,) + 2 



dw 




dxi 





The proof is completed. □ 
Remark 1. By the same way, we can show that, if w £ H^{Q) and r > 7r/(2-v/2) then 

2 

w{xi,X2) cos(aixi) cos{a2X2)dx 

i?2\B{0,r) 

This result improves immediately the results of 



, ^ 16^/27r ,, „2 
da < \\w 



Proof of theorem 2 

Proof. Recall that q, 5e,Re are defined by ([7]), and Gj{Ie), are defined by Lemma 4. For 
aU j G {1, 2}, we define fje : ^ R 

MO = 4^ / G',(/,)(a)e*(«-")da 

Applying Lemma 4, we have Gj{Ie) G ^^(fl^) n L'^{R'^) , so /^^ G C(i?2) n ^^(i?^) and 
Fifje) = Gjile). Applying Lemma 4 again, for all e G (0, e~^), we get 



l-6q 
2 



F(/,,) - F(/,exO < X{B{0, Re))CoRee 

+ x{R\B{0,Re))\F{fje.] 



(17) 
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where Cq is a positive constant depending only I ex- It implies that 



Fifje) - Hfje.) < 2x(i3(0, Re))CiRy-'''l 

+4x(S,) \\f)eAlvn)+'^x{R'\BiO,Re)) F{1 



Ijex, 



Hence, 



Fifje) - FU.ex) < 2CoVi?^ei-6'? + 4m(B,) \\f, 



da 



R'^\B{0,Re) 

It is obvious that IC^irR^e^-^'^ < i?^^ for e > small enough. Moreover, since 

Be C ({a G 5(0, i?,), |L»i(/e.)(a)| < e"} U {a G B(0,i?,), |L>2(/exO(«)l < e^) 

, we apply Lemma 2 (with (/Jo(i) = ^exiT — t)) to get that m{Bi,) < 2R~^ for e > small 
enough. Thus, for e > small enough, we get 



F{f,e)-F{f,e 



Fifjex) 



dadf3 



R^\B{0,R^) 



By Parseval equality, we have 



\fje /jex||/,2(f^) < 



fjs fjex 



2 _ 1 



Fifje) -F{f, ex) 



L2(i?2) 



< 



47r2 



^ + ^ ll/jex|li2(Q) + 2 



V 



F{1 



jex, 



(18) 



da 



ii2\B(0,iJe) 

for e > small enough. Since F{fjex) G L'^{R'^), we obtain that 

limjfje- fjex\\L2^n) =0 
If /jex G H^{Q) then using (jlSp and Lemma 6, we get 

72^/27r 



2 1 / 1 8 2 



||/jex||//i(Q) 



< (ee ||/iex||^i(o) + = 63er (^66 ||/iex||^i(o) + ■ 

for e > small enough. This complete the proof. 



1 ln(ln(e-i)) 
ln(£-i) 



□ 



Proof of Theorem 3 



17 



Proof. We shall use the notations of the proof of Theorem 2. Notice that the assumtion 



R2 



fjex{x).cos{a ■ x)dx 



da < oo, 



is equivalent to F{fjex) G L\R^). Since fjex,F{fjex) G {R^) r\ L'^ {R^) , we get 

JjexiO = ^ / F{J,ex){a)e^^''<Ua 



ij2 

Therefore, 

From (jl7p . we have 



< 



Fifje) - Fifjex) 



Fifje) - F{f,ex) ^^^^^^ < CoTTRte— + 2m(i?,) + 



F{fjex) 



(19) 



do 



For e > small enough, we have CqttR^e 2 ' < i?^ ^ and m{B^) < 2R^ ^. Thus, from (fT9]) . 
for e > small enough, for all j G {1,2}, we get 



1 4 



Since F{fjex) G L^iR'^), we obtain that lirn 



F{fjex) da 

R2\B{0,R^) 

fjexW T,oo(<^) = for all i G {1, 2}. □ 



/je J J ex 1 1 2,00 (Q) 



Remark 2. 14^e can replace R^ defined by (7) 6|/ 

i?, = 10(ln(e-i))^/^° 
to construct a better regularized solution in the case that e is not too small. 

4. A numerical experience 

Assume that T = 1, n= 1/12, A = -1/8. 

We consider the exact data lex = ^jUo,Uq,ut) given by 

f = y sm(7rt), 

TT 

Xi = — sinfvrt). [sin(27rx2)?T-i + sin(47rxi)n2] , 
6 
vr 

X2 = — sinfvrt). [sin(27rxi)n2 + sin (47r2;2) ni] , 
6 

uq = ut = (0,0), 

Uq = (7rsin(47rxi) sin(27rx2),vrsin(27rxi) sin(47rx2)) ■ 



18 



Then the corresponding exact solution of the system (1) — (4) is 

Uex = (sin(7rt) sin(47ra;i) sin(27rx2), sin(7ri) sin(47ra;i) sin(27rx2)) , 
fex = (cos(27ra;i) cos(47ra;2), cos(47ra;i) cos(27ra;2)) • 

For each n = 1, 2, 3, we consider the inexact data In = {(fn, X'''' ,Uo, u*q''iUt) given by 

TT 

X" = Xi+ sin(7rt). [sin(2n7ra;2)rai + 2sin(2ra7ra;i)ra2] , 

TT 

X2 = X2 -\- J— sin(7rt). [sin(2n7rxi)n2 + 2sin(2n7rx2)ni] , 

u^ = uJ^ = (0,0), 

TT 

Un^ = Uq -\ sin(2ra7ra;i) sin(2ra7ra;2) (1, 1) • 

nyn 

Then the corresponding disturbed solution of the system (1) — (4) is 
= Uex H 1= sm(irt) sin(2ra7ra;i) sin(2n7ra;2) (1, 1) , 



fdi ~ fex + 

We get 



(3 3 \ fn 
-y/n ;= I sin(2n7rxi) sin(2n7rx2) H cos(2n7rxi) cos(2n7rx2) 
2 n^^n J 2 



(1,1)- 



</5n = V'; 

\\X^{t,.)-Xf{t,.) 



Uq = Uq, Uj. = UT, 

lho,"--S,IL.(a) = ^'V,E{l,2}, 

and 

II „ _ ||2 _ 5 9,9 

Wljdi JjexWLHn)- 4n 4n3' 
Hence, when n is large, a small error of data will cause a large error of solution. It show 
that the problem is ill— posed, and a regularization is necessary. 

We shall construct the regularized solution as in Theorem 1 corresponding e = n"^/^. 
Prom the straightforward calculation, we obtain that 



32.'sin(JgL)sin(Jg|; 

lap - 247r2^ . ( lal^ - 127r2 



gi{In)i(x) = D{In){a) X (sin(ai)sin(Q!2) - (1 - cos(ai))(l - cos(a2)))x 

2aiQ;2 ^ ^/ri,{aia2 + 127r2(2 - n^)) 



(af _4^2)(ct2 - I67r2) (af - 4n^Tr'^){al - in-^TT^ 
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Thus, the regularized solution defined by 



.n / X 1 f gi{In){oi).D{In){a) 



B{0,Rn) 

where 

For example, if e = 10~^ then 

n = 10^5„ = 0.01389495494, ^„ = 39.52948133, 

and we have some figures about the exact solution fiex, the disturbed solution /"^- and the 
regularized solution/j},g. 



Exact solution 
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Disturbed solution 




Figure 2. The disturbed solution. 



FouriartransFarm dI 9xact solution 




Figure 3. The Fourier transform of the exact solution. 
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Fourier iransform of ragularizad aolulion 




Figure 4. The Fourier transform of the regularized solution. 
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